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Abstract 

Since the invention of space-time coding numerous algebraic methods have been applied to code design. In 
particular algebraic number theory and central simple algebras have been at the forefront of the research. 

In the first part of the paper we will push this direction further and show how the error probability of algebraic 
codes is tied to some central aspects of algebraic number theory and central simple algebras. In particular we prove 
how the error probability of several algebraic codes is tied to the corresponding zeta functions and unit groups. 

In the second part of this paper we turn to study what information theory can say about algebra. We will first 
derive some corollaries from the diversity-multiplexing gain tradeoff (DMT) Zheng and Tse and later show how these 
results can be used to analyze the unit group of orders of certain division algebras. 

Index Terms 

division algebra, space-time block codes (STBCs), multiple-input multiple-output (MIMO), unit group, zeta 
functions, diversity-multiplexing gain trade-off (DMT), algebra, number theory. 

I. Introduction 

In this paper we study the connections between theory of multiple-input multiple-output (MIMO) fading channels 
and algebras, especially algebraic number fields and central simple algebras. We will prove several relations between 
the information theory of MIMO channels and lattice codes that are tighter than in the case of classical Gaussian 
channels. 

The capacity results of the Gaussian channel consider families of codes, where we let the code length grow to 
infinity. As a consequence, we cannot expect that the capacity of the Gaussian channel can tell something non- 
trivial about the geometric structure of a lattice code with finite length, as in such case we are far from capacity. 
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This should not come as a surprise as the capacity proofs use sphere packing arguments. In this sense classical 
information theory, when turned into expressions in lattice theory, is more or less a part of classical lattice theory 
in a weak form. 

However, the theory of fading channels is totally different. When looked as a classical compound channel, the 
quasi-static Rayleigh fading channel has capacity 0. The diversity-multiplexing gain tradeoff (DMT) lH can then be 
seen as a substitute for the classical ergodic capacity results. The difference to the classical case is that the statements 
do consider codes with finite length, and in particular we can consider the DMT of a single lattice space-time code. 
This creates a far tighter relation between lattice theory and information theory of fading channels. This gives a 
good hope that information theoretic results might be turned into non-trivial statements in lattice theory and algebra. 

In this paper we will show how there indeed exists a two way link between multiple antenna information theory 
and lattice theory and algebra. The connecting link is the inverse determinant sums to be presented later. We will 
show that these sums allow us to get a better picture of the error probability of algebraic codes and to apply 
information theoretic results to the study of algebra. 

A. A short description of the content 

We begin by recalling some results of lattices and of DMT. In Section|III]we first show how the inverse determinant 
sums arise naturally from the context of DMT. We then prove some elementary upper and lower bounds for the 
asymptotic behavior of these sums. The upper bound is simply the number of lattice points within a sphere and 
the lower bound comes from the Epstein zeta function. We then derive a stronger lower bound for the asymptotic 
growth of matrix lattices with high dimension. As an warming-up example, we study the determinant sum related 
to the Alamouti code |3| and recognize that it is the truncated Epstein zeta function. Finally, following the union 
bound approach we we give a proof for the fact that the Alamouti code is DMT optimal for a single-antenna 
receiver 

In Section |IV] we begin to study the diagonal MISO codes from algebraic number fields. Here we see how the 
corresponding inverse determinant sum can be asymptoticly approximated by combining the information of the 
geometric structure of the unit group and of the truncated Dedekind zeta function at integer point. As far as we 
know, this is the first time that this connection is pointed out and successfully exploited in the analysis of the error 
probabilities of number field codes. 

In Section [V] we relate the error probability of the division algebra based codes to the Solomon zeta functions 
and to the density of unit group in orders. We then apply the DMT-based lower bound for the inverse determinant 
sums to give density results for the unit group of a family of division algebras. Again the connection with zeta 
functions of orders and unit group is pointed out here for the first time. 

B. Related work 

The study of inverse determinant sums is a natural question in multiple antenna fading channels. For example, 
the work Q by Tavildar and Viswanath is closely related to our approach. However, their work did not consider 
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determinant sums, but eventually restricted their attention to coding schemes, where elementary combinatorial 
methods could be applied. In |]6| the authors studied the blind detection of QAM and PAM symbols. In their 
analysis they applied the Dedekind zeta-function of the field In Example 14.11 we discuss briefly how their 

approach can be seen as the most simple case of our theory. 

The closest (and independent) line of research that is related to our work is done recently by F. Oggier and 
J.-C. Belfiore. In Q they consider Rayleigh fast fading wiretap channels and number field codes. In particular by 
measuring error probabilities in wiretap channel they end up with the same number field sums as we do. In O 
Belfiore and Oggier consider the Rayleigh fading MIMO wiretap channel, where their work also leads to the same 
inverse determinant sums. However, their analysis considers only Alamouti code. 

In the crossroad of ours and the work of Oggier and Belfiore is the work by HoUanti and Viterbo ||9l- They 
considered the error probability of wiretap codes using similar methods to ours. In particular their goal has been 
to improve on the bound given in Section IIV-DI 

C. Main contributions of this paper 

The contributions of this paper are the following. 

• A connection between error probability, Dedekind zeta function and unit group of algebraic number field codes. 

• A formal definition of inverse determinant sums. 

• A connection between Epstein zeta function and inverse determinant sums. 

• General upper and lower bounds for inverse determinant sums. 

• A connection between MIMO error probability, Solomon zeta function and unit groups of orders of division 
algebras. 

• An application of DMT to the analysis of unit groups of orders in Q(z)-central division algebras. 

In Sections |lll] and |IV] we prove the DMT optimality of the Alamouti code and diagonal number field codes, 
when received with a single antenna. Although these results are previously known, the proofs of these results more 
or less rely on recognizing the MISO channel as a parallel channel, and therefore give little insight how the codes 
should be optimized and which code we should use. On the contrary, our union bound approach will offer a good 
insight into the performance of these codes and methods for optimization. 

From the mathematical perspective the question of inverse determinant sums is a special case of general problems 
considering sums over lattice points. However, it appears that general tools, for example from analytic number theory, 
cannot be used effectively here. This is due to the fact that these determinant sums are very sensitive to the choice 
of lattices. 

From another point of view, these sums are related to the Dedekind and Solomon zeta functions. But in our 
case, instead of summing over all different ideals we will sum over all elements, which obviously usually cause 
non-convergence. As far as we know, there is no previous work on the asymptotic behavior of these natural sums. 

In Section [V] we apply information theoretic results to analyze the unit groups of orders of Q(i)-central division 
algebras. This problem can be reduced to a point-counting problem of lattice subgroups in the Lie groups. This is 
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a rather recent but highly developed mathematical area having a rich spectrum of general methods. For the most 
recent approach based on ergodic methods we refer to the monograph by Gorodnik and Nevo ifTOl . The case we are 
considering in Section IV] leads to a point counting of lattice subgroups inside the special linear group SLn{C). The 
case 5i„(R) was already solved in ifTTl . Recently Gorodnik and Weiss lfT2l and Maucourant ifTsl gave very general 
methods for these asymptotic problems, but no explicit formula. Yet, with some algebraic work, their methods will 
give us an exact answer (with small error term). Therefore, it is obvious that our lower bound will not give any 
new information to this problem. However, we think that it is worth pointing out that we can easily say something 
of this highly non-trivial problem, by using information theoretic methods. This suggests that this connection might 
be worth of a further serious study. 

II. The Players 

A. Diversity multiplexing gain trade-off 

Let us now consider a slow fading channel, where we have Ut transmit and receive antennas, and where the 
decoding delay is T time units. The channel equation can be now be written as 

Y = J—HX + N, 
V nt 

where H € M„,,x)it(C) is the channel matrix whose entries are independent identically distributed (i.i.d.) zero- 
mean complex circular symmetric Gaussian random variables with variance 1, and e A/„^xt(C) is the noise 
matrix whose entries are i.i.d. zero-mean complex circular symmetric Gaussian random variables with variance 1. 
X S Mntxxi'C) is the transmitted codeword, and p presents the signal to noise ratio. 

Definition 2.1: A space-time code C is a set of (rit x T) complex matrices. We often use the abbreviation STBC 
for space-time block code. 

Let us suppose we have a coding scheme, where for each value of p we have a code C(p) having |C(p)| matrices 
in M„xt(C). The rate R{p) for code C(p) is then log (|C(p)|)/T in bits per channel use, where the logarithm is 
taken to base 2. Assuming that the scheme fulfills the average power constraint 

I ^'^'^ X<£C{p) 

we then have the following definition from l?). 

Definition 2.2: The scheme C{p) is said to achieve spatial multiplexing gain r and diversity gain d if the data 
rate 

hm - — -- = r 

p->-oo l0g(/9j 



and the average error probability 



lim i^^li^ . ^d. 



p-s-oo log(p) 

In the following we will frequently use the "dotted" notation, where for real valued functions / and g we write 

f{M) = giM), 
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if 



log(/(A/)) _ ^.^ log(ff(M)) 



M->-oo log(M) M-j-oo log(M) 

and similarly > and <. 

Remark 2.1: We will also use f{M) > g{M) in the case where g{M) might not have a limit. By this we only 
mean that g is lower bounded by some function c{]\I) where c(i\/) = f{M). 

B. Matrix Lattices and spherically shaped coding schemes 

In this paper, we will consider space-time codes with rit = T = n, and therefore these codes live in the space 
M„ {€.). Since we like to use results from classical lattice theory in M", without any doubts, the following discussion 
shows how this demand gives us a natural inner product that induces the Frobenius norm in space A/„(C). 

We can flatten X e A/„(C) to obtain a 2n^-dimensional real vector x by first forming a vector of length out 
of the entries (e.g. row by row, or vectorizing that is column by column) and then by replacing each complex entry 
with the pair formed by its real and imaginary parts. This defines a mapping a from A/„(C) to R^"^: 

a: X ^ x = a{X) (2) 

which is clearly R-linear: 

a{rX + r'X') = ra{X) + r'a{X'), r, r' S M. (3) 



Let ||X||^ = y^Ti{XJX) denote the Frobenius norm of X. Note that the following equality holds: 

11X11 



J2J2\^rs\'-\HX)\\E, (4) 

i=l j = l 



where \\-\\^ denotes the Euclidean norm of a vector. This makes a an isometry. It also gives us a natural inner 
product in the space M„(C). Let us suppose we have matrices X,Y E A/„(C). We can then define {X,Y) = 
^{Tr{XY^) = {a{X),a{Y)), where the last () is the natural Euclidean inner product in R2"^ 
Definition 2.3: A space-time lattice code C C A/„(C) has the form 

where the matrices Bi, . . . ,Bk are linearly independent over R, i.e., form a lattice basis, and k is called the rank 
or the dimension of the lattice. 

Definition 2.4: 1241 If the minimum determinant of the lattice L C Ai„(C) is non-zero, i.e. it satisfies 

inf |det(X)|>0, 

we say that the code has a non-vanishing determinant (NVD). 

We now consider a coding scheme based on a fc-dimensional lattice L inside i\/„(C). For a given positive real 
number M the finite code is 

L{M) = {a\a(^ L,\\a\\p < A/}. 
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In the following we will also use notation 

B{M) ^ {a \ a e MniC) \\a\\p < M}, 

for the sphere with radius M. 

The following result is well known, but can also be proved by using Proposition lA.ll in Appendix with s = 0. 

Lemma 2.1 (Spherical shaping): Let L be a /c-dimensional lattice in Af„(C) and L{M) be defined as above; 
then 

\L{M)\ = cM^ + o{M^-^), 

where c is some real constant, independent of M. 

Proof: This follows directly from Proposition lA.ll ■ 
In particular, it follows that we can choose real constants Ki and A'2 such that 

KiM^ > \L{M)\ > K2M''. 

The coding scheme, where the finite codes are sets 

= (5) 

yields the desired number of codewords for multiplexing gain r and for each p level. Furthermore, the sets Cl [p) 
clearly do satisfy the average energy constraint ([T|). Henceforth we simply ignore the scaling factor of SNR in 
the channel equation as it is irrelevant in DMT calculations. 

Remark 2.2: One can wonder if the peak power energy normalization (|5]l will lead to too low average energy 
for DMT considerations. It is indeed true that this methods does result into lower codeword energy than bounding 
the average energy as in ([T]l. However, we can show that this difference is irrelevant in the DMT context. Let us 
suppose we have a fc-dimensional lattice in A/„(C) and measure the average energy of the code p^^L (p^) 

According to Proposition I A. 1 1 we have 

On the other hand we also have that \L{p^)\ ^ p''". Combining these, we can see that from the DMT perspective 
it is irrelevant which normaUzation we use. 

III. Inverse Determinant Sums Over Matrix Lattices 

In this section we consider the inverse determinant sums over matrix lattices and show how these sums create a 
connection between algebra and information theory. We first begin with a non-rigorous introduction, which shows 
how these sums appear naturally as a continuation of more familiar sums. 

Consider a /c-dimensional lattice code L{M) C C" for the additive complex Gaussian noise channel 

y — x + n 
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where x_ £ L{M) and n is a length-n complex Gaussian random vector with zero mean and covariance matrix /„. 

We have the famiUar expression of the pairwise-wise error probabihty (PEP) upper bound for confusing x_ to 
in the receiver 

If the codewords from the code L{M) are sent equiprobably, we can upper bound the average error probability by 
the following sum 

sGi,0<||^||^<2A/ 

where the term 2M follows as we have to consider differences of codewords. The right-hand-side is indeed a well 
known exponential sum taking values on lattice points. 

Let us now consider a quasi static Rayleigh fading channel with single transmit and rir receive antennas. Assuming 
the channel matrix is known perfectly to the receiver but not to the transmitter, in this case we have for the code 

L{M) e C" 

P{x^^)<- 

\\x-x\\e 

and the corresponding sum is 



x'^L,0<\\x\\j^<2M \\ — \\e 

We can then see that if 2nr > k the right-hand-side is the truncated Epstein zeta function lll4l . 

We now turn to the more general case of having a fc-dimensional NVD lattice L{M) C A/„(C) for a slow 
Rayleigh fading MIMO channel with n transmit and n,. receive antennas 

Y = HX + N, 

where H and W are respectively the channel and noise matrices and where X E L{M). In terms of PEP, we have 
for X :=^X' 

P(X X') < 

and the corresponding sum is 

< V i . 

- ^ |det(X)|2». 

XGL,0<||X||y<2M ' ^ ^' 

We summarize the three cases above below. 

• Single antenna channel AWGN; P^. is upper bounded by the sum of e""-"^^ an exponential sum. 

• Single antenna slow fading channel: Pg is upper bounded by the sum of .i , an Epstein zeta function. 

• Quasi-Static Rayleigh fading MIMO channel: Pg is upper bounded by the sum of | dct(X)P'"- ' ™ inverse 
determinant sum. 

We will see that the behavior of this last sum is the most peculiar. While in the first and second cases we can make 
the sums converge by increasing iir, whilst in the last case of inverse determinant sums we will show that it might 
not converge. 
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A. Basic Problem 

Let us suppose that we have a fc-dimensional lattice L C A/„(C). The finite codes with multiplexing gain r 
derived from the spherical coding scheme are 

Assume there are receive antennas. Then following the union bound together with PEP based determinant 
inequality ifTSl . we get the following bound for average error probability for code Cl{p^"^'') 

P < V ^ ^ — 

^ , |dct(X)|2". ' 
XGL(2p™/'') ' ^ 

where in the above sum (and also in the rest of the paper) the all-zero matrix is always excluded from the sum. To 
simplify the notation, we will be considering the following sum 

1 



XeL{M) ' ^ " 

B. Elementary and information theoretic bounds for inverse determinant sums 

In this section we provide some general bounds that, give us some ideas of general behavior of inverse determinant 
sums. They can be applied to the performance analysis of codes as well as to algebra, as suggested. First we give 
rather elementary bounds in Proposition 13.31 and then in Proposition 13.41 we apply DMT results of Zheng and Tse 
to get a far stronger lower bound. This is in strong contrast to the results of Proposition 13.31 

Lemma 3.1: [16| Let A and B be invertible matrices in A/„(C) and let ai > • • • > a„ be the eigenvalues of 
AA^ and &i < • • • < be the eigenvalues of BB^ . We then have that 

n 

2 



\AB\\l > ^a,6,. 



4=1 

Lemma 3.2: Let us suppose that X is a set of matrices in A/„(C) and that A is an invertible matrix in A/„(C). 
If / is such a function that 

\B{M) r\x\< f{M), yM 

then there is such a constant K that 

\B{M)nAX\ < f{KM),yM, 

where AX = {AX : X e X}. 

Proof: Let us suppose that Ai is the smallest eigenvalue of A^A. According to Lemma [TTI we now have that 
for all the elements AX e AX, \\AX\\p > Ai \\X\fp. It follows that for matrix AX, where 

UXy < M, 

we must have that < -^=- We now see that is a suitable constant for K. 
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Armed with the above lemma, we are now ready to provide some easy upper and lower bounds for the asymptotic 
behavior of the sum X]xgl(a/) | dct(x)|"' ^'^^ ^ being a fc-dimensional NVD-lattice L in A/„(C). 
Proposition 3.3: Let us suppose that L is a /c-dimensional NVD-lattice in Af„(C), with 



We then have that 



mindet(L) := min |dct(X)| = 1. 
o^xeL 



KM'' > y , ^ }^^, > A^iM'^-™", k~mn>0 

XeL{M) ' ^ 

KM'' > V ^——>K2\og{M),k-mn = Q 

KM'' > y , ^ > ifg, k-mn<0, 

- ^ dct(X - 



for some constants A', K\, Ki, and A'3. 

Proof: Minkowski inequality gives us that 

\X\ 



|det(X)| < 
We then have that 



F 



> — mn 

1 . sr^ \/n 



^ idct(x)|™ - „ „ ^ ii^iir' 

X£L(M) ' ^ ^' \\X\\p<R,X£L " 

By applying Proposition lA. 1 1 we get the lower bounds. 



On the other hand, let us now consider the worst case and suppose that | det(X)| = 1 for all nonzero X ^ L. 
In this case we have 

y I , = y 1 = |L(A/)| < iVA/^ 

^ dct(X) ™ ^ I V ;i - 

where iV is a constant independent of M and where the last inequality follows from Lemma 12.11 

■ 

We next focus the situation where L is a 27T,^-dimensional lattice in Ai^r, (C), i.e., k = 2v? . It turns out that the 
DMT results by Zheng and Tse can be used to provide a far stronger lower bound on the inverse determinant sum 

Proposition 3.4: Let L be a 2n^ -dimensional NVD-lattice in A^„(C) with. We then have that 

XeL{AI) ' ^ 

Proof: Consider the previously mentioned coding scheme Cl{p) = L (^p^^ with k = 2n^. Just as 
previously, the union bound gives us that 

1 



, |det(X)|2". ■ 

The optimal DMT given by Zheng and Tse, on the other hand, gives us that for integer values of r we have that 
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It follows that 

for integer values of r. Let us now set R = 2p''/^" and we then have that 

The maximum here is achieved obviously with r = 0, but in this case we do not have growth for our matrix sum. 
The next integer point is r = 1. In this case we have 

SUR) > 

■ 

We now see that for 2n^-dimensional lattices the asymptotic behavior of Sl{M) is around M^" (or worse). 
The most interesting thing here is that no matter how large rir we choose this result is valid. We also see that in 
some sense the behavior of the sum is almost the worst possible. 

Remark 3.1: As far as we know there is no existing Uterature for these natural sums. 

C. The inverse determinant sum and DMT of Alamouti code 

In this section we will show that Alamouti code does reach the bound in Proposition 13. 3 1 This result then allows 
us to rediscover the DMT of Alamouti code when received with rir antennas. 
Let us use the following notation 

' xi + X2i -{x^ + x^i)* 
.xs + Xii {xi + X2i)* 
The corresponding lattice of Alamouti code can be written as 

Lalarn = 0, 0, 0) + • • • + ZA{0, 0, 1, 0) + ZA{0, 0, 0, 1), 

which is a 4-dimensional lattice in 7\/2(C). We then consider the corresponding inverse determinant sum 

1 



A{xi,X2,X3,X4) 



^ |det(X)|2'"- 



Proposition 3.5: We have that 



E ^i^r(:^^^^i«g(M), 

where Ki and K2 are some constants. 

Proof: Due to the orthogonality of the rows of the Alamouti code we have that for any codeword in Laiam 



|det(X)| 

We now have that 



V2 



E \r\r,i-lY\\2m E 



■)2m 
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The rest follows from Proposition lA.ll ■ 
Remark 3.2: In particular if m is large enough the inverse determinant sum of the Alamouti code is the Epstein 
zeta function. 

Corollary 3.6: When received with rir antennas, the Alamouti code achieves the optimal DMT curve 

(r, 2nr{l-r)+) 

for any lattice codes with 4-dimensions in M2(C). 

Proof: In order to study the DMT of the lattice code Laiam, we use the previously introduced coding scheme 
Laiamip^^'^)p^^^^ ■ The usual union bound argument now gives us the following bound for the error probability of 
making a mistake in reception when transmitting an arbitrary codeword 



According to Proposition 13.61 we then have that 

Pe < p-2"'-(i^'')^(21og(2p'-/2)) 

where A' is some constant independent of p. This gives us that Alamouti code does achieve the claimed DMT. 
Proposition 3.3 in IT] shows that this upper bound is the best achievable for lattice codes with 4-dimensions in 
A/2(C). ■ 

IV. Algebraic number field codes, Dedekind zeta function, density of the unit group, and DMT 

Let us now consider inverse determinant sums arising from algebraic number field codes iflTl . In particular we 
will show how the error probability of these codes is tied to the unit group and Dedekind zeta function of the 
corresponding algebraic number field. These connections allow us to give a better look at the behavior of these 
codes and to prove their DMT optimality. For simplicity let us consider a degree n cyclic number field extension 
K/Q{i), where the Galois group is (cr) = {ai, ■ ■ ■ ,ct„}, and Ok is principal ideal domain (PID). See Section 
IIV-DI for comments about these conditions. 

We can define a relative canonical embedding of K into Af„(C) by 

■0(2;) = diag(cri(x), . . . ,CT„(a;)), 
where x is an element in K. The ring of algebraic integers Ok has a Z-basis W = {wi, . . . , W2n] and therefore 

^{Ok) = yj{wi)Z + ■■■ + V(W2„)Z, 

is a 2n-dimensional lattice of matrices in Af„(C). The main reason to use such construction is that for each element 
nonzero a G Ok, we have that | dct{ip{a))\ > 1. Let us now suppose that we have a 2n-dimensional number field 
lattice code L C A/„(C) and that we are considering the coding scheme, where the finite codes are chosen by the 
method of Lemma ISTI Before measuring the DMT for this type of codes, we will need some concepts and lemmas. 
We begin with an example. 
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Example 4.1: The simplest example of the previous construction arises from the trivial extension Q(i)/Q(i). 
The Galois group then consists simply of the identity element. We then have a lattice L = Z[i] £ C, which is a 
2-dimensional lattice in A/i(C). The corresponding sum is 

y y ' 



I 2nr- 



We can now see that this is actually the truncated Epstein zeta function and we are again reaching the bound in 
Proposition 13.31 However, we can look at this problem also from another angle that will be easier to generalize. 
We notice that for every element in a; G Z[i], we have |xp = ^Q(i)/Q(2;). 

We know that Z[i] has only 4 invertible elements {1, —1, i, — i} and that Z[i] is a principal ideal domain. Therefore 
for every ideal a;Z[i], we have exactly 4 different generators {x, —x, ix, —ix}. We can the write 

y y 4 

x6L(M),x#0 11-^11 N{I)<VM 

which is the truncated Dedekind zeta function CQ(i)('S) at point s = n,.. In particular when we let grow to 
infinity we get that the sum Y.xeL(M),x^o approaches 4CQ(,)(nr). 

We point out that this approach was earlier taken in |j6l. Yet one can use it only in the case when the extension 
has degree 1. We will next show how this can be extended to more general number fields. 

Consider a cyclic extension K/Q{i), where [K : Q{i)] = n. The corresponding inverse determinant-sum 



U{x)\\p<M,x^Ok ' \\i,{x)\\p<M,x<^OK ' '^^ x£X(M) ' '^'^^ 

where X{M) is a set of elements x E Ok, \\''P{x)\\p < M , each generating a separate ideal in Ok- Accordingly, 
Ax{M) is the number of elements x £ Ok, ||V''(2;)IIf — ^ ' ^^^^ generating the same principal ideal in Ok- 

If we neglect for the moment the terms Ax{M), and consider only sum J^xex iNK/vix)]"'- ' ^^^^ ^ 

part of Dedekind zeta function (^k at point n^. 

In the following we will give bounds for Ax{M). The bounds depend only on the value of M; they are independent 
of the choices of x. 

A. A bound for Ax{M) 

We will begin our analysis with Ai{M). This is simply the part of the unit group of the ring Ok, consisting 
of elements u e O}^ that ip{u) G B{M) C A/„(C). 

Lemma 4.1: Let K/Q(i) be a cyclic field extension with [K : Q{i)] = n. Then the set 

Ai{M) = {^{u)\u e 0*K, mu)\\p < M}, 

has an upper bound 

AiiM) \AiiM)\ < N{logMf-\ 
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where is a constant independent of M. 

Proof: For ease of reading, the proof to this lemma is relegated to Appendix |B] ■ 
Based on Lemma 1431 we can upper bound the value of Ax{M) for all x. 

Proposition 4.2: Let /i/Q(z) be a cyclic field extension with \K : = n and let x E Ok be a non-zero 

element with ||'!/'(2:)IIf — Then 

A.,{M) = \{u\U(xu)\\p<M,ueO*K}\<N'{\ogMr-\ 

where N' is a constant independent of M as well as of the element x. 

Proof: Given x E Ok, we can write tp{x) ~ diag(xi, . . . , x„). The condition ||-!/'(a:;)|lp < M implies \xi\ < M 
for all i. We also have that [xij • • • \xn\ > 1- It follows that for all i 



kil > (6) 



Now for any u e O^^- is a unit such that j|7/'(Ma;)||^ = \\tp(u)tlj{x)\\ ^ = ||diag(a;iui, . . . , x„u„)j|^ < A/, where 
Ui = (Ti{u). Since ux E Ok, we have that |a;i||ui| < M for all i and (|6]l implies \ui\ < M". Therefore we have 
that < y/nAr\ Lemma 143] now gives that 

A^{M) < Ai (V^A/") < N (log(v^Ar))""^ < N' (log Af)""^ , 

where N' is a constant independent of M. ■ 
The essential part of this result is that we could find a constant K, such that A'(logA/)"^^ upper-bounds all 

the A^{M), for all x G Ok, with \\i}{x)\\p < M. 

B. A bound for partial zeta sum 

In the following we will denote by I^- the set of integral ideals of the ring Ok- 
Proposition 4.3: Let K/Q{i) be a cyclic field extension with [K : Q{i)] — n. Then 

xex{M) ' 

where X{M) is a subset of Ok in which each element x generates a separate integral ideal and satisfies ||V'(2;)llir ^ 
M, as defined in dlVb . is a constant independent of M 

Proof: Using basic properties of algebraic norm and AM-GM inequality we have 

\det{^{xW = \NK/q{x)\ <mx)\\l\ 

for any element x € X{M) C Ok- This implies 

\Nk/q(x)\<M''" 

where X is a set of elements x G Ok, each generating a separate ideal in Ok 
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( E ^] <(log(M^"))^ 



From the relation between ideals and element norms we can further upper bound the above quantity by 

|7Vk/q(x)|<A/2" |JVk/q(/)|<A/^" 

where / represents an integral ideal. Note that the right-hand-side corresponds exactly to the beginning of the 
Dedekind zeta function at point Ur- It then follows that 

y I < ' 

|JVjf/Q(/)|<A/"" 

where the first inequality is based on a similar reasoning as in IfTSl Prop. 7.2, Cor 3] as well as some elementary 
approximation. ■ 
We remark that if > 1 the bound in Proposition 14.31 is tighter See Subsection IIV-DI for a discussion. 

C. The inverse determinant sum and DMT of algebraic number field codes 

Armed with Lemma 14.21 and Proposition 14.31 we are now ready to continue the derivation of ( IIVI ) to obtain an 
upper bound for the DMT of number field codes. 

Proposition 4.4: Let K/Q{i) be a cyclic field extension with [K : Q{i)] = n. Then 



where N is some constant independent of M. 
Proof: 

Continuing from (lIVI i we have 



< Ni [logMf-^ N2 (logM)'" , 

where the first inequality follows from Lemma l4~2l to upper-bound Ax{M) with a constant A^i, and the second 
inequality is due to Proposition 14.31 with another constant ■ 
Remark 4.1: Here we have an example of a 2n-dimensional lattices, where the growth of inverse determinant 
sums is logarithmic. Comparing this bound to that in Proposition 13.31 we can see that we are somewhat close to 
lower bounds if rij. = 1, but are not even in correct class when > 2. This suggests that the bound in Proposition 
3.3| is not very tight. 

Finally, we are ready to determine the DMT curve achieved by the previously discussed number field code L, 
by which we mean the following. Let K/Q{i) be a cyclic field extension of degree n; then 

L = {diag(cri(x), • • • ,CT„(a;)) : xgOk} 

which is a 277-dimensional lattice. Given SNR p and multiplexing gain r, let 

CLip^) = p-iL{pi). 
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be the corresponding finite code obtained by the spherical coding scheme 

Theorem 4.5: If the receiver has antennas, the code Cl{p^) achieves the following DMT curve 

(r, — r)^). 

Proof: Note that L is an NVD lattice. It can be easily shown that the maximal pair-wise error probability 
achieved by Cl{p^) is == p-^^^ii-r) hence Pg > p-"">-(i-i-) y p^j. ^j^g uppgj- bound on Pe, the usual union 
bound argument gives 

„—nr.n(l — r} — n,.Ti(l — r) 

^- ^ , |det(X)|2»- - |det(^(a;)|2". U^gl^P ij -P 

xeL{2p2) \[4,{x)\\p<2p2 ,xeOK 



where the last inequality follows from Proposition 14.41 Combining the upper and lower bounds on proves the 
claim. ■ 

D. A remark on the constant values 

In Proposition 14.41 we showed the following result 

U{x)\\p<M,xeOK ' 

For cyclic extensions, where Ok is PID, we point out the assumption of being cyclic is not needed anywhere. This 
bound is also true in the case where Ok is not PID, but it is only looser 

Yet, this general result can (and will) later be applied also to the division algebra case, but due to the generality 
lacks tightness. Below Let us shortly describe how our methods can give quite tight asymptotic bounds for number 
field codes, when the number of receiving antennas is greater than 1. 

We note that the term (log7\/)^" coming from Proposition 14.31 can simply be replaced with C-ffC^-r), when 
is greater than 1. This already reduces the bound in Proposition 14.41 to A^(log A/)"^^. We can say furthermore a 
few words about the constant N . 

Main theorem in |!20'| gives us the following bound that is true for large enough M 

R(n — Ij! 

where lo is the number of roots of unity in K, R is the regulator of the number field K and > 1 is an arbitrary 
constant. A short study of the proof of Lemma 14.11 gives us that actually 

A,{M)<A^{M). 

Assume > 1. By collecting all the previous results we now have 

U{x)\\^<M.x<^Ok ' 



j , iv± la loi^v.^ v^iiwugii aiiu U/ ^ J- an 

probably quite tight asymptotically for the leading term (log(M))""^), but generally we are overestimating because 



where N = "^"„_2)! — > -^^ is large enough and d > 1 is an arbitrary constant. For the PID case this bound is 



by using Dedekind zeta function we have included in the sum all the ideal classes that might not be principal. 
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V. An Application of Information Theory of Fading channels to Algebra 



In the previous section we connected Dedekind zeta function and the density of a unit group of an algebraic 
number field to the error probability of algebraic codes. In this section we will prove that if we consider space 
time-codes from orders of Q(i)-central division algebras we will get an analogous result. Here the ring of algebraic 
integers Ok will be replaced with an order A and Dedekind zeta function with Solomon zeta function. We would 
also need an analogous result for the density of unit group of an order of a division algebra. 

However, in the case of division algebras the corresponding density result is highly non-trivial. It can be seen 
as a question of arithmetic lattices in Lie groups. Explicit formulas do exist for the case of special linear group 
S'L„(R) ifTTl . In our case the relevant Lie group is S'L„(C) and there do exist general methods lfT3l . lfT2l . but no 
explicit formulas. Therefore, instead of using algebraic results to reveal something about information theory, we 
will investigate what information theory can say about algebra. 

We will show how Proposition 13 .41 can be used to derive some asymptotic results of the density of units of orders 
inside division algebras. Due to the nature of our methods we cannot expect that we could achieve better results 
than those by using general methods from |fT3l and lfT2ll . However, the achieved result do suggest that Proposition 
3.41 carries some non-trivial information of the behavior of inverse determinant sums and that there might exist a 
more direct link connecting outage probability and Lie groups. 



A. Division algebras. Orders, Zeta functions and Unit Groups 

Let E/Q{i) be a cyclic extension of degree n with with Galois group G{E/Q{i)) ~ (cr). Define a cyclic algebra 

V = {E/Q{i), a,-f) = EOuE® u^E © • • • © u^'^E, 

where u G 2? is an auxiliary generating element subject to the relations xu = u(j{x) for a\\ x £ E and u" = 7 G 
(Q(i))*. We assume that 2? is a division algebra. 

Definition 5.1: A Z-order A in P is a subring of 2?, having the same identity element as T), and such that A is 
a finitely generated module over Z and generates 2? as a linear space over Q. 

We next show how an order can be turned into a matrix lattice with NVD-property. Consider 2? as a right vector 



,Tl— 1 



space over E and every element x = xo+ uxi 

ip{x) 

Xo 7(T(a:„_i) "fa'^{xn-2) ■■■ "f(T'^~^{xi)^ 

Xl cr(xo) 7(T^(x„-i) 7o-""^(a;2) 



Xn-i <= 2? has the following representation as a matrix 



/ 



X2 



cr(xi) 



a'^ixo) 



"fa" 



\Xn-i a-{x„-2) a'^^Xn-a) 
A simple and easily describable order is the natural order 



(9) 



A„ 



)uOe®u'^Oe< 



17 



where is the ring of algebraic integers in E. 

From our perspective the most important properties of these Z-orders are the following. If A is a Z-order in a 
division algebra V, then tp{A) is 2n^-dimensional NVD lattice in 7\jf„(C), with 

\dct{X)\ > 1, 

for all the nonzero elements X in ip{A). 

We need a few more concepts. Let V be an index-n Q(i)-central division algebra and A be a Z order in T>. The 
(left) Solomon zeta function l23l of the order A is 

Ca(s) ^ ^ 



leiA 

where 5fts > 1 and Ia is the set of left ideals of A. The fact we need from this function is that it is indeed a 



converging series 

The result connecting this sum to our matrix lattice consideration is the following 

|det(?/>(x))|2" = [A : Ax]. (10) 

The unit group A* of an order A consists of elements a; e A such that there exists an y e A with xy = 1. The 
unit group O'^ of the ring of algebraic integers Ob is very well known and has a simple structure. However, this 
is not the case for the group A*. In most cases it is extremely mystical Il22ll . 



B. Inverse determinant sums and density result on A* 

As already seen, the lattice structure of V'(A) allows us to consider sum 

which follows from an consideration of error probability of lattice codes in fading channels. 

In this section the main result is Proposition |5.3l but we first need some preliminaries. First, we have the following 
result, which is an analogy with Proposition 14.21 in orders. 

Lemma 5.1: Let V be an index-n Q(i)-central division algebra and A be a Z order in V. For any x E A with 

0< <M, 

A.,{M) = \{u\ mxu)\\p <M,ueA*}\< \ij{A*) n i?(Af")|. 

Proof: Let Ai, . . . , A„ be the eigenvalues of ip{x)ip{xy . The condition ||-(/'(a;)|| < AI implies < AP for 
all i. We also have that |Ai| • • • |A„| > 1 by the NVD property. It now follows that 

I^^I^^IT^rryV. (11) 

For any unit u E A* such that |l?/'(itx)||^ = \\t/j{u)ilj{x)\\ p < M, Lemma [TT] shows 

\\t{j{u)il}{x)\\p > min{Ai, • • • , A„} • 
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Combining (fTTT l and ||V;(M)-0(a;)||^ < M gives < M", and hence proves the claim. ■ 

The following result is analogous with Proposition 14.41 The only thing preventing us for giving a numerical 

bound is the missing information about the unit group, i.e. the value of |V'(A*) H i?(Af")|. 

Proposition 5.2: Let V be an index n Q(i)-central division algebra and A be a Z order in V. Then, for Ur > 2n 

E |detW..))P». ^^'-^(^^)^^(^^")'' 

Uix)\\j,<M.xeA ' 



where is a constant independent of M. 
Proof: The sum 



can be written as 

Ax{M) 



E 



where X{M) is some collection of elements x G A, ||'0(2;)IIf — ^^'-h generating a separate (left) ideal. 
Accordingly, the number Ax{M) is the number of elements inside •ip{A) n B{M) that generate the same ideal Ax. 
Set rir — sn with s > 2 by assumption. We then have 

^^^^^^|det(^(x))p- ^^A. [A : A.]- 

is a part of the Solomon zeta function of order A at point s > 2. Therefore, the above sum is always bounded by 
some constant N independent of M. 

From the ideal theory of orders we have that if Aa; = Ky, then x and y must differ by a unit. Therefore, Lemma 
Ocan be applied to show Ax{M) < |^(A*) n B(M")| for all x G A with < ||V'(a;)||;^ < M. Finally, we have 
the following series of inequalities 

^ |det(^(x))p- - [A-A^^ iaTJ^] <N\^iA)nBiM )\, 

and the proposition is established. ■ 
We now combine this result with Proposition 13.41 

Proposition 5.3: Let A be an order in an index n = 2m Q(z)-central division algebra V. Then, 

|V'(A*)nB(M)|>M2"-2_ 

Proof: Note that ijj{A) is a 2n2-dimensional lattice in A/„(C). By Proposition 13.41 

V >m2"^-2" (12) 

^ |det(V'(x))|2".- ^ 



holds for any > 1. On the other hand. Proposition 15.21 gives 



^ det(V' X P"- - 1*^^ ^ ^ ''I' 

x<^A,mx)\\p<M I ''^ 
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for rir > 2n and for some constant N independent of M. Combining the above two results gives 

|V'(A*)nS(M)|>M2"-2_ 

■ 

This simply means that for big enough values of M the hypersphere B{M) with radius M in M„(C) contains 
at least close to Af^"~^ elements of ^'(A*). On the other hand, ^/'(A) has approximately M^"^ elements inside the 
same hypersphere by Lemma 12.11 While the number of units is small compared to the whole number of points of 
the lattice, it is still remarkably larger than in the case of number fields where it is in class (logAf)""^. 

Remark 5.1: We remark that (fT2l i does already tell that an order must include a "plenty" of elements with small 
determinant value. This is a direct consequence of Proposition 13.41 It is somewhat surprising that such a general 
result can give us some understanding of this very special problem. 

VI. Discussion and further study 

It is clear that to get a better picture of the performance of space-time codes, we should be interested not only 
in the minimum determinant of the lattice, but also about the spread of determinants in the lattice. The asymptotic 
study of inverse determinant sums does provide us with a well-defined tool to get a rough picture of the spread 
of determinants of space-time lattice codes. Proposition 13 . 3 1 can then be seen as a benchmark, when studying such 
sums. The analysis performed for number field and division algebra codes does prove that it is indeed possible to 
estimate these sums for practical codes. We hope that the tools and connections we developed in algebraic sections 
will be beneficial to the future study of algebraic space-time coding. 

Proposition 13.41 is a relatively direct consequence of the DMT by Zheng and Tse. In particular 13.41 follows 
directly from the asymptotic outage bound in ||4| . The connection between the growth of inverse determinant sums 
and outage bound does follow from union bound and therefore is probably deemed to be loose. However, it seems 
to be able to tell something non-trivial about the spread of norms of elements in the orders of division algebras 
and in particular about the unit groups. 

There are several new directions this research can lead to. Let us describe two most obvious ones. 

It is quite obvious that our general results on the asymptotic behavior of inverse determinant sums are not very 
strong yet. The bounds we gave are based on a very crude approximation. An interesting question is the following. 
If Proposition 13.41 does provide some information of the behavior of determinant sums that is not achievable with 
more direct methods and if so, can this approach be simplified and generalized. This study could also lead to a 
new code design criteria that could maybe replace DMT, at least for lattice codes. 

The other, more practical issue is to analyze the inverse determinant sums of algebraic codes more thoroughly. 
Here the first question is to measure the correct growth class, polynomial, logarithmic, etc. It is also crucial to 
determine the constants. Of course, before the analysis, we have to normalize the fundamental parallelotopes of the 
corresponding lattices. We hope that this approach will bring us a surprise or two. 
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Appendix A 
A USEFUL Result for Summing over Lattices 

It is hard to find a exact reference for the following well known result. Because it is so central in DMT context 
(see Lemma 12.11 Remark 12.21 and Proposition 13.31 ). we will sketch a proof for it. 
Proposition A.l: Let L be a fc-dimensional lattice in W\ Then 

HiM''+'' < \\x\\% < i^2^/"+^ s+k>o 

xeL{M) 

i^3log(^/)< ME<Hi\og{M), s + k^Q 

xeL{M) 

J2 M%<H4, s + k<0, 

xeL(M) 

where Hi are constants independent of M. 

Proof: The first step is to move from a fc-dimensional lattice in R" to a fc-dimensional lattice in M'''. This can be 
done because instead of considering lattice L, we can just as well consider lattice AL C M'', where A e A/fexri(R) 
is any matrix such that AL is a fc-dimensional lattice in M'^'. Let A,ni„ and A,nax be respectively the minimal and 
maximal eigenvalue of Al^ A, which are obviously independent of M. Clearly, Amax > Amin > 0. The following 
inequaUty is obvious for any x & L 

min 

It then follows that for any s > (for s < we get similar inequalities) 



This implies that if we can prove the results for the fc-dimensional integer lattice inside R'^, then the claims are 
true for any fc-dimensional lattice in R". Specifically, we can choose A so that it will take L to Z'^ c R'" with the 
set of basis vectors 

{(1,0...,0),(0,1,...,0),...,(0,0,...,1)}CR^ 

Therefore, instead of an arbitrary fc-dimensional lattice in R" we can just as well consider lattice Z'' C R'^. 

Let 1 = 1^ C W". When s < -fc the sum is simply a truncated version of the Epstein zeta function and it is 
known that the sum converges. Let us now suppose that s > — fc. The following inequalities follow from repeated 
application of Euler-McLaurin summation formula ( 1251 p. 17) 



Jl<\\y\\E<M xeL{M) 

<K2f miEdy, 

Jl<\\y\\^<M+l 

where Ki are constants independent of M. If we suppose that s 7^ — fc, we have 



\\y\\Edv~y2M 



fc+s 
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and the result follows. 

Let us now suppose that s = — fc. We then have that 

which gives us the second result. 

■ 

Appendix B 
Proof pf Lemma |4~T] 

Lemma 14.11 is an elementary corollary to Dirichlet unit theorem, but we give a proof, as it sheds some light on 
the question. 

Proof of Lemma \4~l\ The number field K has signature (0,n). The Dirichtlet unit theorem tells us that the 
unit group O}^ has the following multiplicative structure 

= Ufree X Uroots = Z X Urootsj 

where Uroots is a finite torsion group containing roots of unity in Ok- Let us consider the mapping / : M" 

f{u) = (log|CTi(u)|,log|cr2(w)|,...,log|CT„(w)|). 

It is well known that f {Ufree) is a (n — 1) -dimensional lattice inside R". 

Let us now consider il;{Ufree) H B{M). If i^lu) happens to be inside a ball B(M) with radius AI we have in 
particular that |cri(M)| < AI for all i. It follows that for coordinates ai{u) having absolute value greater than 1 we 
have log(|cri(w)|) < log(A/). On the other hand if |cr(iii)| < 1, we have that |/og(|(Ti(u)|)| < (n— l)log{AI), which 
is a consequence of the facts that for positive coordinates log(|t7i(u)|) < log(Af) and Yll=i log(l'''i(^)l) — 0- 

Therefore, if '4}{u) is inside a ball of radius AI then f{u) is inside a cube with side of length {n — 1) log(Af). 
We have that f {Ufree) is a (n — 1) -dimensional lattice, and therefore a cube with {n — l)log(A/) side has less 
than A^log(A/)"^^ elements, where iV is a constant independent of AI. It follows that \4'{Ufree) H B{M)\ < 
N \og{AI)"^^ . Now each of the elements in ip{Uroots) is a unitary matrix. Hence for any u ~ UfUr G O}^ with 
Uf £ Ufree and Ur G Uroots, wc havc \\il;{ufUr)\\p = \\ip{uf)\\p. Therefore we see that 

\^{Ufree X Uroots) H B{AI)\ - \^{Ufree) H B{AI)\ ■ \Uroots\- 

It follows that 

mUfree X Uroots) B{AI)\ < N \og{AI)^-^ ■ \Uroots\. 

As the group Uroots is finite the claim follows. 
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